Number of observable features in the acoustic-Raman spectra of nanocrystals 
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Low-frequency Raman scattering spectra are presented for gold nanocrystals with diameters 3.5 
and 13 nm. The frequencies of the Raman peaks but also their number are shown to vary with the 
nanocrystal size. These results are analyzed using both the continuous elastic medium approximation 
and an atomistic approach. We show that the number of atoms in the nanocrystal determines 
an upper limit of the number of observable Raman features. The frequency range in which the 
continuous elastic medium approximation is valid is defined by comparison with the calculations 
based on the atomistic approach. 
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Many advances have been made in the last decades re- 
garding the experimental investigations and theoretical 
descriptions of confined acoustic vibrations in nanocrys- 
tals (NCs). In particular, many works have been devoted 
to inelastic light scattering spectra (Raman or Brillouin) 
which are known to depend on the size, shape, compo- 
sition, environment and lattice structure of the NCs^i 2 - 
Recent papers have reported a variety of spectra for very 
small (a few nanometers) to quite large (hundreds of 
nanometers) NCs, especially for metallic NCs. For NCs 
sizes smaller than the optical wavelength, the Raman se- 
lection rules dramatically restrict the number of Raman 
active vibration modest For larger nanoparticles (hun- 
dreds of nanometers), several vibration modes can con- 
tribute to the Raman/Brillouin spectral However, while 
the assignment of the spectral features has received a lot 
of attention, the question of their number remains largely 
unaddressed. 

Assuming a continuous and isotropic medium, the vi- 
bration eigenmodes of a spherical particle can be derived 
using Lamb's model.— One obtains spheroidal and tor- 
sional modes here labeled Sf and T™, respectively; I and 
n being the angular momentum and overtone index. The 
number of vibration eigenmodes is infinite. This predic- 
tion is however bound to fail for clusters consisting of few 
atoms. 

In this work, we present experimental results obtained 
for gold NCs which demonstrate that the number of 
acoustic vibration modes and associated Raman peaks 
strongly depend on the NC size, i.e., on the number of 
atoms in the NC. Assuming a continuous elastic medium, 
we derive an upper limit for the number of observable 
Raman peaks. In order to properly describe the high 
frequency vibration modes with wavelengths compara- 
ble to the interatomic distance, we also use an atomistic 
description. In this approach, the non-linearity of the 
acoustic phonon dispersions as well as the contribution 
of surface atoms to the vibrational dynamics are taken 
into account. 



A. Raman setup, samples and spectra 

Low-frequency Raman scattering was measured on ap- 
proximately spherical 13 and 3.5 nm gold NCs, chem- 
ically synthesized in colloidal solutions. In all Raman 
experiments the Argon/Krypton ion laser intensity was 
around 0.5 mW//zm 2 at the focal plane. The scattered 
light was collected through a 50 x objective of a confocal 
microscope with a 0.8 numerical aperture and dispersed 
using a Jobin-Yvon T64000 spectrometer which allows 
for high rejection of the Rayleigh scattering. The laser 
spot size is about 10 /^m 2 and the incidence angle is 60° 
with respect to the normal of the sample surface. 

In order to enhance the Raman scattering from the 
13 nm NCs, the latter were coated on gold nanodisks 
and the solvent was evaporated. The NCs were then 
excited at 647.1 nm in resonance with the surface plas- 
mons of the nanodisks as described in Ref. 0. The nan- 
odisks provide a spectrally sharp surface plasmon res- 
onance associated with spatially localized electric near- 
fields which allows for a strong enhancement of the Ra- 
man scattering. As shown in Fig. Q] several overtones 
of the NCs vibration modes are observed^ Three addi- 
tional low-frequency peaks are present but not resolved in 
this spectrum due to the Rayleigh scattering. They have 
been reported elsewhere and assigned to the S\ (split into 
Eg + T2 g by elastic anisotropy) and Sf vibrations^ 

The Raman spectrum of the 3.5 nm NCs is also shown 
in Fig. [TJ The NCs were coated on a 300 nm Si0 2 /Si 
substrate and excited at 520 nm in resonance with both 
the surface plasmons of the nanoparticles and one of the 
Si02/Si optical cavity modes in order to take advantage 
of a double resonance effect: surface plasmons and inter- 
ference enhanced Raman scattering^ This is particularly 
useful in the case of small particles with a very weak op- 
tical density. Because of the smaller size, the lowest fre- 
quency peaks are now observed (S% split into E g + Ti g 
modes) together with the S\ band. This modes assign- 
ment is based on the Raman selection rules (breathing 
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FIG. 1. (Color online) Low-frequency Raman spectra of gold 
NCs with diameters of 3.5 (top) and 13 nm (bottom). The 
Raman shift v has been multiplied by the diameter d and 
divided by the average longitudinal sound velocity vl- The 
arrows indicate the spectral features due to S\ (large arrows, 
red online) with E g and T2g and Sq vibrations (small arrows, 
black online) with n = 1, ... 5 from left to right. 

















♦ 


♦ ♦ 















12345678 




12345678 



diameter (nm) 

FIG. 2. Top: variation of the highest frequency of the acoustic 
vibrations of a gold NC as a function of its diameter. Bottom: 
maximum number of breathing vibrations Sq for a gold NC 
as a function of its diameter. 



modes generate only polarized Raman scattering whereas 
quadrupolar modes give rise to both polarized and depo- 
larized scattering). 

In Fig. [TJ the Raman shift has been multiplied by the 
diameter d of the NCs in order to compare the number 
of spectral features for different NC sizes. In addition, it 
has been normalized to the longitudinal sound velocity 
of gold. In that way, the frequencies of the breathing 
vibration modes (calculated using the continuous elas- 
tic medium approximation) are located approximately 
at integer abscissa (short arrows in Fig. [1]). This figure 
clearly shows that changing the size of the NCs signifi- 
cantly changes the number of Raman features: the Ra- 
man spectrum of the 13 nm NCs exhibits more spectral 
features than that of the 3.5 nm NCs. 



B. Lamb's model 

We first discuss the vibrational dynamics of spheri- 
cal NCs in the frame of Lamb's model. 5 In this model, 
the vibrating medium is assumed to be continuous and 
isotropic. The vibration cigcnmodcs arc obtained from 
the Navier-Stokes hydrodynamic equations using Vl — 
3330 and vt — 1250 m/s for the longitudinal and trans- 
verse sound speeds of gold;2 According to the Raman se- 
lection rules^ for NC sizes small compared to the optical 
wavelength, only spheroidal vibration modes with angu- 
lar intergers £ = and 1 = 2 (i.e., Sq and S% modes) are 
Raman active in agreement with experiments^ In this 
work, we will focus on the existence of breathing vibra- 
tion modes Sq and related Raman features. 

It is well-known 10 that the continuous elastic medium 
approximation becomes questionable when the NC con- 
tains only few atoms and/or the wavelength of the vi- 
bration mode is comparable to the interatomic distance. 



In order to take into account these limitations, a sim- 
ple approach consists in considering only the 3./V lowest- 
frequency vibration modes (N being the number of atoms 
in the NC) with angular periodicity larger than the in- 
teratomic distance (£ max = tt— ^ — This simple 
rule gives a rough estimation of the number Nb rea thing of 
allowed vibration modes and their frequency range v max . 

Fig. [2] shows Nbreathing and v max as a function of the 
NC diameter. It appears that a gold NC can sustain an 
Sq breathing mode only if the NC diameter is larger than 
1.0, 1.9, 2.8, 3.7 and 4.6 nm for n = 1, 2, ... 5 respectively. 
Therefore, the number of observed Raman peaks in Fig[T] 
is determined (at least partly) by the NC size. For in- 
stance, the Sq vibration mode is observed in the Raman 
spectrum of the 13 nm NC (Fig. [T]). But such a mode is 
absent in the spectrum of 3.5 nm NCs because it does not 
exist for this NC size (independently from the associated 
Raman scattering efficiency). 

Unlike Nbreathing, v max weakly depends on the NC size 
and is between 100 and 125 cm -1 for d ranging from 1 
to 8 nm. For large diameters, this value is close to the 
Debye frequency for gold (118 cm -1 ) due to the similar 
assumptions used in Lamb and Debye models. Hence, the 
frequency range where the low-frequency Raman peaks 
are expected weakly depends on the NC size. It certainly 
cannot reach very high frequencies for small sizes 



C. Atomistic approach 

The previous discussion provides a basic approach to 
the understanding of the number of vibration modes a 
NC can sustain. However, since this approach is based 
on the assumption of a continuous and isotropic elastic 
medium, it fails to correctly describe the vibration dy- 
namics of surface atomsi^ and the non-linearity of the 
phonon dispersions. As a matter of fact, in bulk gold the 
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FIG. 3. (Color online) Projection of the overtones of th 
breathing vibration modes Sq with n = 1 to 5 from top t 
bottom on the atomistic eigenmodes of a gold NC with 75 
atoms (d ~ 3.16 nm) as a function of the atomistic eigenfrs 
quencies. The frequencies of the Sq vibrations are indicate 
by vertical lines (red online) where n — 1, ... 4 from left t 
right. 



phonon dispersion significantly deviates from the linea 
variation around 100 and 60 cm" 1 for longitudinal an 
transverse phonons, respectively 

Atomistic calculations have been carried out usin 
the same method and parameters as those reported i: 
Ref. [TJ. The size imposed limitation of allowed vibrs 
tion modes was already pointed out for the Sq modei- 
In the following, we shall focus on the overtones Sq c 
the breathing mode for the NC size d ~ 3.5 nm inves 
tigated experimentally in this work. Calculations wer 
performed for gold NCs with 755, 959 and 1270 atom 
corresponding to d ~ 3.16, 3.34 and 3.7 nm, respective!} . 
We obtain v max ~ 150 cm -1 a value close to the one 
estimated using Lamb's model (Fig. [2]). 

Following Ref. Q, the atomistic eigenmodes have been 
projected onto Lamb modes Sq in order to evidence some 
of the atomistic vibrations which may contribute to the 
Raman spectra. Figs. |31 IH and [5] show such projections 
for NC with 755, 959 and 1270 atoms. Due to the or- 
thogonality and normalization of the atomistic modes, 
the sum of the squared projections onto a given Lamb 
mode is equal to one. 

The Sq and Sq projection peaks are very close to the 
frequencies of the Sq and Sq Lamb modes. This points 
out the good agreement between the vibration frequen- 
cies calculated using the atomistic model and Lamb's 
model. The frequency of the Sq projection peak is signif- 
icantly lower than that of the Sq Lamb mode, mainly 
because of the non-linearity of the phonon dispersion 
not accounted for in the continuous medium approxima- 
tion. In addition, significant projections are observed in 
a broad frequency range, especially for the smallest NC 
size (755 atoms, Fig. [3]). Significant projections occur 
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FIG. 5. (Color online) Same as Fig.EJfor a gold NC with 1270 
atoms (d ~ 3.7 nm). 

also at the frequencies of the Sq and Lamb modes. 
This is even more pronounced for the projections onto 
the So and Sq Lamb modes (Figs. EH HI andO, which 
indicates that the atomistic eigenmodes strongly differ 
from the Sq and Sq Lamb modes. 

Fig. IH1 shows the vibrational density of states (VDOS) 
calculated using the atomistic model and Lamb model. 
For the latter, the size-imposed restriction of the number 
of allowed vibration modes has been taken into account 
as well as the 2£+ 1 degeneracy of the spheroidal and tor- 
sional modes. The shape of the atomistic VDOS is similar 
to the one already reported for bulk gold^ whereas the 
Lamb VDOS is quite different. First, it is not smooth 
especially at high frequency due to the large degeneracy 
of the involved modes. Second there is no high frequency 
tail but rather an abrupt cut-off (due to the size-limited 
number of modes). It clearly appears that Lamb model 
is not reliable in the high frequency range of the acoustic 
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FIG. 6. (Color online) VDOS for a 3.34 nm (959 atoms) Au 
NC calculated from Lamb's model (left) and the atomistic 
model (right). The bin size is 3 cm -1 . 

spectrum. 

In order to investigate more precisely the validity of 
Lamb model, we plot in Fig. [7] the frequencies of the n- 
breathing-like atomistic modes as a function of a pseudo- 
wave vector Q n . We define the n-breathing-like atomistic 
modes as the modes with projections on Lamb modes Sq 
greater than 75% of the maximum of these projections 
(see Figs. 131 IH and [5]). To these modes, we associate a 

pseudo-wave vector Q n = — j^- 2 - that accounts for the 
radial dependence of their displacement field; vs™ being 
the frequency of the Sq mode and vi the longitudinal 
sound speed in gold. In Lamb's theory, the displacement 
field of the Sq Lamb mode is proportional to ji(Q n r )l ji 
being the spherical Bessel function of order 1.— 

Fig. [7] shows the frequencies of the n-breathing-like 
atomistic modes as a function of Q n . They are com- 
pared to the linear dependence of the Lamb frequencies 
i/gy. Despite a relative dispersion of the atomistic fre- 
quencies (especially for Q n > 5 nm -1 ) due to the uneven 
projections, the frequencies of the atomistic modes ob- 
tained for the different NC sizes approximately overlap. 
As can be seen, for Q n < 5 nm , the vibration frequen- 
cies obtained using the atomistic and Lamb models are 
in a very good agreement. However, for Q n > 5 nm -1 , 
Lamb model clearly overestimates the vibration frequen- 
cies in comparison with the atomistic model. Since Q n 
is a rough estimation of a pseudo-wave vector for the 
spherical symmetry, the plots in Fig. [7] reflect the phonon 
dispersion curves for both models and therefore point 
out the difference between the atomistic calculations and 
Lamb model. 



D. Comparison with experiments 

1. Raman scattering 

The spectra shown in Fig. [T] clearly demonstrate the 
main result of our calculations, namely that the number 
of observable Raman peaks is determined by the NC size. 
However, one must first check that the inhomogeneous 
broadening due to the size distribution of the NCs does 
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FIG. 7. (Color online) Frequencies of the n-breathing-like 
atomistic (red) and Lamb 5ff (black) modes as a function of 
Q„. The number of atoms in the Au NC is 755, 959 and 1270 
for the circles, triangles and crosses respectively. 



not prevent the observation of some Raman features. 

Assuming that the linewidth of the Raman lines is 
due only to inhomegeneous broadening (i.e., neglecting 
the homogeneous contribution) and that the scattered 
intensity for each peak does not depend on size, the 
linewidth/frequency ratio is the same for all the Raman 
peaks. This ratio is r ~ 0.35 for the S\ E g peak of the 
3.5 nm NCs and r ~ 0.25 for the Sq peak of the 13 nm 
NCs. Two consecutive overtones will significantly overlap 
when their frequency separation becomes equal to their 
full width at half maximum. Because ~ nwoi, this 

occurs for n 1/r i.e., n > 3 and 4 for the 3.5 and 13 
nm NCs, respectively. If the homogeneous broadening is 
not negligible (which is the case in this work since the 
width of the size distribution is about 15% of the aver- 
age size for both samples) or if the scattered intensity 
for a given peak depends strongly on the size then the 
linewidth increases more slowly with n and more over- 
tones can be observed. Therefore the size distribution 
does not explain the non-observation of the Sq peak in 
the spectrum of the smallest NCs. However it may be 
responsible for the non-observation of the Sq peak of the 
13 nm NCs. 

The Raman scattering efficiency decreases with in- 
creasing overtones^ which makes the observation of such 
high order modes (n > 1) very challenging,- The Ra- 
man spectrum of the 3.5 nm NCs (Fig [lj does not ev- 
idence any peak corresponding to the first harmonic of 
the breathing mode (Sq) despite the resonance and in- 
terference enhancement. However the atomistic (Figs. [3l 
IU [5]) and Lamb (Fig [2]) models point to the existence of 
such a mode. Modeling of the Raman scattering efficien- 
cies based on an atomistic description of the NC vibra- 
tion eigenmodes and their interactions with the electronic 
states is required in order to more precisely address this 
issued 
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2. Transient absorption 

The breathing modes of gold NCs can also be observed 
using time-resolved transient absorption experiments. 
Up to now, fewer overtones of the breathing vibrations 
have been reported using this technique compared to 
Raman scatteringi&iZ Recently, Juve et al*£ conducted 
time-resolved transient absorption experiments on small 
(1.3 < d < 3 nm ) platinum NCs and investigated the 
size-dependence of the fundamental Sq breathing modes: 
no discrepancy between the measured vibration frequen- 
cies and the frequencies calculated using the Lamb model 
has been found for NCs as small as 1.3 nm. 

According to our atomistic calculations (Fig. [7]), a 
strong deviation of Lamb frequencies from the atomistic 
frequencies is expected for Q n > 5 nm -1 . Assuming that 
the elastic properties of platinium and gold are similar, 
Q n > 5 nm -1 corresponds, for the fundamental breath- 
ing mode n = 1, to a NC diameter d < 1.2 nm. Below 
this critical size, Lamb model fails to describe the vibra- 
tion dynamics of the NCs. The minimum size (1.3 nm) in- 
vestigated by Juve et a/.— is close but still above the criti- 
cal size, thus explaining the good agreement between the 
measured frequencies and those estimated using Lamb 
model. 

Time-resolved transient absorption experiments^ per- 
formed on very small Au NCs (1.1 < d < 2.2 nm) have 
evidenced an acoustic vibration mode, apparently the Sq 
breathing mode, whose frequency does not change with 
the NC size. This behavior is consistent with our atom- 
istic model and in particular with Fig. [71 The largest NC 
size (d = 2.2 nm) for which this frequency independent 
mode is observed^ is in quite good agreement with our 
predictions. However, the frequency (around 70 cm -1 ) 
measured for the smallest NCs is smaller than the one ex- 
pected from the atomistic model. Note that this value is 
close to the first peak in the atomistic VDOS as pointed 
in Ref. [TH, suggesting that the observed mode may not be 
the Sq breathing mode. This again points to the need for 
an improved model based on an atomistic description of 
the NC and an accurate model of the electron-vibration 
coupling. 

E. Other nanoparticles 

While this work was focused on gold NCs, similar re- 
sults are expected for other metallic and semiconductor 



NCs. However, a special feature of gold is the large lon- 
gitudinal/transverse sound velocity ratio. The frequency 
difference between two consecutive breathing modes is 
proportional to vt, while this difference for other vibra- 
tions symmetry (S^o anc ^ ^7) 1S mahily proportional to 
Vt- Hence, there are more modes having a lower fre- 
quency than an Sq vibration mode in gold than in most 
of the other materials. For example, the number of such 
modes is about two times larger for gold than for silver 
and about three time larger than for copper. Hence the 
impact of the size- induced reduction of the number of Sq 
vibration eigenmodes is more easily evidenced in the case 
of gold NCs. 



III. CONCLUSION 

The number of peaks in the inelastic light scattering 
spectra of small gold NCs is shown to depend on the NC 
size. This result can be qualitatively explained by re- 
stricting the number of vibration eigenmodes according 
to simple rules based on the number of atoms in the NC. 
This effect comes in addition to the well-known 1/d de- 
pendence of the acoustic vibration frequencies. Atomistic 
calculations provide a more exact but also more complex 
picture than Lamb's model. Complexity arises from the 
non-linearity of the phonon dispersion curves which is 
very pronounced for vibration modes with wavelengths 
comparable to the interatomic distance. By comparing 
the sequence of frequencies obtained from the atomistic 
and Lamb's models, we were able to define a frequency 
above which the continuous medium approximation is no 
longer valid. This frequency is neither size- nor mode- 
dependent. It is mainly determined by the non-linearity 
of the phonon dispersions. 
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